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Abstract. Two-dimensional colloidal suspensions subject to laser interference patterns with decagonal 
symmetry can form an Archimedean-like tiling phase where rows of squares and triangles order aperiodically 
along one direction [J. Mikhael et aL, Nature 454, 501 (2008)]. In experiments as well as in Monte- 
Carlo and Brownian dynamics simulations, we identify a similar phase when the laser field possesses 
tetradecagonal symmetry We characterize the structure of both Archimedean- like tilings in detail and point 
out how the tilings differ from each other. Furthermore, we also estimate specific particle densities where 
the Archimedean-like tiling phases occur. Finally, using Brownian dynamics simulations we demonstrate 
how phasonic distortions of the decagonal laser field influence the Archimedean-like tiling. In particular, 
the domain size of the tiling can be enlarged by phasonic drifts and constant gradients in the phasonic 
displacement. We demonstrate that the latter occurs when the interfering laser beams are not adjusted 
properly. 



PACS. 61.44.Br Quasicrystals - 82.70.Dd Colloids 



1 Introduction 

Unlike ordinary crystals, quasicrystals possess point-sym- 
metries, e.g., with five- or ten- fold rotational axes, which 
are not allowed in periodic crystals pL|i2j. Nevertheless, 
quasicrystals exhibit Bragg-reflections due to their long- 
range positional and orientational order. Since they show 
exceptional material properties [3 , there is much interest 
in understanding and controlling the growth of quasicrys- 
tals. Therefore, a lot of research has been devoted to study 
the ordering of atoms on quasicrystalline surfaces [4,5 ,6] 
;[ail[9l[IOl[IIl[Il[I3|ll|l5]. However, in these experiments 
. it is very difficult to determine the exact positions of the 
' atoms. Therefore, a model system consisting of micron- 
sized colloidal particles subject to a laser interference pat- 
tern has been used to investigate the dynamics and order- 
ing of particles in a two-dimensional quasicrystalline po- 
tential p!6lll7H 18 ] ll9 l '2Q". Interesting structures have been 
observed, such as phases with 20 bond directions [19 or 
colloidal orderings consisting of rows of squares and rows 
of triangles termed Archimedean-like tilings p!8l [T9] (see 
also Fig.[T]). The latter will be the main topic of this work. 

Interference patterns of five or seven laser beams have 
decagonal or tetradecagonal symmetry, respectively (see 
Fig. [2j). In such laser fields, the colloids are forced to- 
wards the regions with highest laser intensity, i.e., for the 




Fig. 1. Archimedean- like tiling phases induced by (a,b) the 
decagonal and (c,d) the tetradecagonal substrate. (a,c) show 
experimental results, (b,d) are obtained by Monte-Carlo simu- 
lations. The figures show Delaunay triangulations of colloidal 
point patterns. Bonds longer than l.las {l.2as for experi- 
ments) are omitted, where as is the spacing of the parti- 
cles when placed on an ideal triangular lattice. The potential 
strength is Vo/{kBT) = 13 in (a), Vo/iksT) = 10 in (c), and 
Vo/(/cbT) = 20 in the simulations (b,d). The density is given 
by as/av = 0.57 in (a), as/ay — 0.58 in (b), and as/ay = 0.7 
in (c,d). 
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Fig. 2. The external potentials are created by interfering laser beams with a symmetric arrangement. An interference pattern 
with 10- fold rotational symmetry is obtained with five laser beams (a-d), 14- fold symmetry is realized with seven beams (e-h). 
(a,e) show how the laser beams are arranged and (b,f) the wave- vectors projected onto the sample plane. (c,g) demonstrate the 
interference patterns observed in experiments, and in (d,h) they are calculated from Eq. ([T]). The bar in the upper right corners 
of (d,h) denote the length scale av = Stt/ \Gj\, where Gj are the projected wave vectors as shown in (b,f). 



colloidal particles the interference patterns act like a sub- 
strate with quasicrystalline symmetry. If the influence of 
the laser field is weak, the colloidal particles form a trian- 
gular lattice in case of sufficiently strong repulsion forces 
or they are in a fluid phase when colloidal interactions 
are weak. For large laser intensities, the ordering of the 
colloids is determined by the interference pattern, i.e., a 
quasicrystalline phase is observed (see also the phase di- 
agrams in Fig. [3] and in [19 ). In this work we are es- 
pecially interested in phases that occur at intermediate 
laser strengths where there is a competition between the 
repulsive colloidal interactions (that alone would lead to 
periodic triangular ordering) and the external forces due 
to the laser field which favor aperiodic quasicrystalline 
structures. For certain particle densities that we will de- 
termine later in this article, the competition leads to a 
colloidal ordering that consists of rows of triangles and 
rows of squares that establish an aperiodic order in one 
spatial direction, (see Fig. [T]). The structure observed is 
close to an Archimedean tiling, where the rows of squares 
and triangles form a two-dimensional periodic order. It 
was therefore named Archimedean-like tiling when it was 
first realized in experiments using laser fields with decago- 
nal symmetry [18] . Later Archimedean-like tilings also oc- 
cured in simulations that studied the complete phase be- 
havior on a decagonal substrate ([19 , see also Fig. [3]). 



The purpose of this article is to characterize the struc- 
ture of the Archimedean-like tiling phases in detail, to 
discuss the differences of the ordering in laser fields with 
decagonal and tetradecagonal symmetry, and finally to de- 
termine the densities where these phases occur. In addi- 
tion, we study the influence of phasonic distortions on 
colloidal ordering in laser fields with decagonal symmetry. 
Phasons are excitations that only exist in quasicrystals. 
Like phonons, they are hydrodynamic modes, i.e., they do 
not cost free energy in the long- wave- length limit [21^ 22 . 
Phasons are an ongoing main topic in current research on 
quasicrystals and intensively discussed in the field [231I2H 
[25]. 

In Sec. [2] we introduce our system. The details of the 
experimental realization as well as the numerical simu- 
lations are explained in Sec. [3l In Sec. H) we analyze the 
structure of the Archimedean-like tilings and calculate the 
densities where these phases occur. We study the conse- 
quences of phasonic distortions on the Archimedean-like 
tilings in Sec. [5] and finally conclude in Sec. [6l 



2 Model 

A charged stabilized suspension is confined between two 
glass plates and subjected to laser interference patterns 
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potential strength Vo/(kgT) 

Fig. 3. Phase diagram obtained with Monte-Carlo simulations 
for colloids on the decagonal substrate (the complete phase di- 
agrams are presented and explained in [19 ). Unlike quasicrys- 
talline phases that we observed at all densities, Archimedean- 
like tiling structures are only found for intermediate potential 
strengths within a small range of densities (marked green). 



with quasicrystalline symmetries. Five or seven laser beams 
with identical linear polarizations are employed to obtain 
interference patterns with decagonal or tetradecagonal ro- 
tational symmetry, respectively (see Fig.[2j). For vanishing 
tilt angle ^ of the beams, the potential in the xy plane 

is [miis]: 

AT-lAT-l 

^(^) = E E [(G. -Gk)-r + 4>, - 4>kl (1) 

where N is the number of beams, (j)j are the phases of 
the laser-light waves, and Gj their wave- vectors projected 
onto the xy plane [cf. Fig. [J^b) and (f)]. The prefactor is 
chosen such that — Vb gives the minimum value of the po- 
tential. We usually set = for all j corresponding to a 
laser field without any phononic or phasonic displacements 
and a center of exact 10- or 14-fold rotational symmetry 
at r = 0. In addition, for the decagonal potential, we ex- 
plicitly analyze the effect of phasonic distortions on the 
Archimedean-like tiling in the last section. In this case, 
the phases are used to specify the phononic displace- 
ment field u(r, t) = [ux{T^t)^Uy{r^t)] and the phasonic 
field w(r, t) = [wx{r,t)^Wy{T,t)] following the convention 
ofRef. [2ll[22], 

(r, t) = u (r, t) • G^- + w (r, t) ■ Gsj mod 5- (2) 

In section [5] we will show how a uniform phasonic dis- 
placement w(r, t), which increases linearly in time, or how 
a static w(r, t) with a constant gradient can be used to 
align the rows of the Archimedean- like tiling phases along 
specific directions of the underlying substrate potential. 



3.1 Experiment 

The colloidal system we used is an aqueous suspension of 
highly charged sulphate-terminated polystyrene particles 
with a radius of R = 1.45 jam from Interfacial Dynam- 
ics Corporation with an average polydispersity below 4%. 
The pair interaction of particles can be described by a 
screened Coulomb potential as introduced below in Eq. 
(|3]). As sample cell we used a silica glass cuvette with 
200 jam spacing. Before starting the measurement, it was 
connected to a deionization circuit to reduce the salt con- 
centration of the suspension. This circuit also included an 
electrical conductivity probe (typical ionic conductivities 
below 0.08 /iS/cm), a vessel of ion exchange resin, and a 
peristaltic pump. After deionization, the suspension was 
inserted into the cell which was then sealed. Before ap- 
plying the quasiperiodic potential, the density of the par- 
ticles was increased in the field of view. This is achieved 
through thermophoretic and optical forces using vertically 
incident fiber coupled infrared laser with a wavelength of 
A = 1070 nm and maximum output power Pmax = 5 W 
and an argon ion laser (A = 514 nm, Pmax = 7 W). A 
circular monolayer was formed on the lower surface of the 
sample cell with a diameter of about 500 /im. The verti- 
cal fluctuations of the particles are also suppressed by the 
argon ion beam, therefore the system can be considered 
as two-dimensional. The monolayer spontaneously adapts 
the hexagonal packing. After reaching the desired particle 
density, the infrared laser was turned off. The high par- 
ticle density was kept constant using an optical tweezer 
(A = 488 nm, Pmax = 7 W) scanned around the cen- 
tral region of the system, i.e., forming a circular corral. 
The quasiperiodic potential was created by interference 
of five or seven linearly polarized beams of a Nd:YV04 
laser (A = 532 nm, Pmax = 18 W). Due to optical gra- 
dient forces, this interference pattern acts as an external 
potential on the particles. The depth Vq of the singular 
deepest potential well (see Sec. [2]) scales linearly with the 
laser power. Vb was determined by calibrating the laser 
intensity for three laser beams, i.e., when they create a 
triangular potential. The length scale ay was tuned by 
changing the angle of incidence of the laser beams with 
respect to the sample plane. The different colloidal phases 
were identified with the help of a Delaunay triangulation. 
It creates nearest-neighbor bonds between vertices that 
were defined by the maxima of the colloidal density dis- 
tribution. For Archimedean-like tiling structures the bond 
length is a bimodal distribution with the ratio of the peak 
positions close to ^/2. Bonds longer than 1.2a5' were re- 
moved from the triangulation which resulted in the char- 
acteristic square-triangle tiling. 



3 Methods 

In the following, we shortly introduce our methods, which 
we employed to analyze the results from experiments and 
simulations. 



3.2 Simulations 

The interaction between the colloids is modeled by the 
pair potential 

[Z*eexp(/i:i?)]^ exp{-i^d) 

^^'^^"4Wr(l + /^i^)^ d 
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of the Derjaguin-Landau-Verwey-Overbeek theory [271 128], 
where d is the distance between two colloids and the in- 
verse Debye screening length. The prefactor also depends 
on the radius of a colloid, its effective surface charge 
Z*, and the dielectric constant of water e^. 

We quantify the particle density by the spacing as 
of the particles when placed on an ideal triangular lat- 
tice. The occurrence of Archimedean-like tilings mainly 
depends on the density of the system. All other parameters 
can be varied over huge ranges. The main limitation to the 
choice of parameters is that there has to be a competition 
between the colloidal interactions and the interaction with 
the substrate. Therefore, we usually choose the parameters 
such that the density used in the simulations is slightly 
above the density where the triangular to liquid phase 
transition occurs. For example, the parameter set for the 
decagonal substrate employed already in [19] is suitable 
to observe Archimedean-like tilings, i.e., R = 1.2 /im, 
Z* = 1000, Cr = 78, T = 300 K, ay = 5.0 /im, and 
= 0.25 jam. Whereas a parameter set closer to the 
experimental values R= 1.45 /im, Z* = 400000, = 78, 
T — 300 K, ay = 6.5 /im, and = 0.2 /im (cf. [20 ) was 
employed to simulate colloidal ordering in the tetradecago- 
nal laser field. We usually performed Monte-Carlo simula- 
tions using the Metropolis algorithm [30] . Periodic bound- 
ary conditions were implemented for simulation boxes that 
contained about 500 to 1400 colloids. For comparison or to 
study dynamical properties, we also employed Brownian 
dynamics simulations as specified in [17 . 

When domains with Archimedean- like tiling structures 
occur, their rows are aligned along the respective five or 
seven equivalent directions of the laser potential. Usu- 
ally within one simulation box, one finds many domains 
with different orientations. There is a tendency that small 
domains anneal out in time, i.e., adjacent domains grow 
slowly within the simulation process. The structures shown 
in Fig. [1] are about the largest domains we achieved by 
pure Monte- Carlo simulations. However, we discovered dif- 
ferent methods to support the formation of larger do- 
mains. While in Monte- Carlo simulations the direction 
of a domain may change quite easily, such domain flips 
are rare in Brownian dynamics simulations. Therefore, 
Brownian dynamics simulations can help to grow large 
uniform domains. It is also possible to predefine a pre- 
ferred direction. For example, by using small simulation 
boxes where the potential is discontinuous at the bound- 
aries, the rows of an Archimedean-like tiling start to form 
along the boundaries. Furthermore, phasonic distortions 
support the growth of Archimedean-like tilings along spe- 
cific directions as discussed in Sec. [5) 

4 Archimedean-like tilings 

The structures we observe in the two-dimensional qua- 
sicrystalline potentials at certain particle densities consist 
of rows of squares and rows of triangles and therefore are 
similar to the Archimedean tiling shown in Figure SJ In 
general, Archimedean tilings are formed by regular poly- 
gons and are characterized by the number of edges of the 




Fig. 4. Archimedean tiling of type (3^ -4^). 



polygons that meet in each vertex [31 . The tiling shown 
in Figure m is called (3^ • 4^)-Archimedean tiling, because 
going around a vertex, one first finds three regular tri- 
angles and then two squares (see red arrow in Fig. |4]). 
Very recently similar Archimedean tilings were also ob- 
served for binary mixtures of nanoparticles situated be- 
tween a cubic and a dodecagonal phase [32 . Furthermore, 
an Archimedean tiling of type (3^ • 4 • 3 • 4) occured in simu- 
lations of monoatomic systems on periodic substrates [33] 
El[35j[36]. 



4.1 On decagonal substrates 

The Archimedean-like tiling on a decagonal substrate dif- 
fers from the perfect (3^ • 4^)- Archimedean tiling by the 
frequent appearance of double rows of triangles [see Figs. 
[Tl^a,b)]. From the simulation data we were able to identify 
the sequence over a range of —Gay < ^ < lOOay that cor- 
respond to 77 single or double rows of triangles. Within 
this range the observed sequence is exactly the sequence 
of the Fibonacci chain, if one identifies a double row with 
L and a single row with S [cf. Fig. mb)]. A Fibonacci 
chain is obtained by starting with element L and repeat- 
edly applying the iteration rules L LS and S ^ L. 
For example, the first iterations are L, LS^ LSL^ LSLLS^ 
LSLLSLSL, etc.. 

Note, the squares and triangles of the Archimedean- 
like tiling are not perfect regular polygons. A row of tri- 
angles has almost the same height as a row of squares. In 
the simulations, we find the ratio of the mean row height 
ly divided by the mean particle spacing Ix in the direction 
along the rows to be ly/lx ~ 0.95, i.e., the triangles are 
stretched in the direction perpendicular to the rows while 
the squares are stretched in the direction along the rows. 
We will calculate the ratio ly/lx in subsection 14.31 

At a first glance, the Archimedean-like tiling structure 
seems to be periodic along the rows of triangles. How- 
ever, some aperiodic modulations along the rows can be 
observed. For example, in the structures shown in Figs. 
ma,b) only very few bonds are perfectly horizontal. Whereas 
most of the bonds between two neighboring rows of trian- 
gles are almost horizontal, the direction of all other bonds 
deviate from the horizontal direction in an aperiodic way. 
Furthermore, there is a modulation of the bond lengths 
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Fig. 5. Colloidal order close to symmetry centers in simu- 
lations: (a) On decagonal substrates double rows of triangles 
occur at symmetry centers, (b) in tetradecagonal laser fields 
double rows of squares are found. 



along the direction of the rows. As a result, some rows 
of squares exist along a direction that is rotated by 27r/5 
from the main horizontal row direction. 



4.2 On tetradecagonal substrates 

The Archimedean-like tiling phase on a tetradecagonal 
substrate consists of large periodic regions that correspond 
to a perfect Archimedean tiling of type (3^ -4^). The 
periodicity is only interrupted by a few double rows of 
squares. Fig. [5jb) shows the position of the colloids on 
the substrate. Obviously, a double row of squares occurs 
at a local symmetry center of the laser field. This bright 
spot corresponds to a deep potential minimum that is sur- 
rounded by 14 shallower potential wells. It therefore con- 
stitutes the center of a region with 14-fold rotational sym- 
metry. Since in a tetradecagonal laser field local symme- 
try centers are rare [20 , only few double rows of squares 
exist in the Archimedean-like tiling structure. Fig. [5fa) 
demonstrates that for the decagonal substrate double rows 
of triangles occur at symmetry centers. Since there are 
much more symmetry centers in a decagonal laser field 
[20] , much more double rows exist in the Archimedean- 
like tiling phase of a decagonal substrate compared to a 
tetradecagonal laser field. 



4.3 Densities with Archimedean-like tilings 

In the following, we want to estimate the densities where 
Archimedean-like tiling structures occur. At a suitable 
density, the rows of colloids should be oriented along lines 
with many minima of the potential landscape. In the decago- 
nal laser field, we can determine the distances between 
such lines by averaging the potential V^(r) along the re- 
direction, i.e., 



25 



-5-2 cos {27ry/Xi) - 2 cos (27r?//A2)] , (4) 



with Ai 



and A2 





Fig. 6. We calculate the mean bond length in horizontal direc- 
tion Ix in a way that many colloids lie along the 27r/5-direction 
(circled red). The mean row hight ly is chosen to fit one of the 
two length scales of the substrate in ^-direction. 



We now calculate the density of an Archimedean-like tiling 
whose mean row height ly fits to one of the lengths Ai or 

A2. 

The simulation results show that the mean bond length 
Ix in the direction of the rows differs from the row height 
ly. Therefore, in the following, we aim at calculating the 
ratio ly/lx' So far we have determined the row height ly 
of an Archimedean- like tiling that fits to the decagonal 
substrate. To estimate the bond length Ix^ we argue as 
follows. The rows are oriented along a symmetry direction 
of the decagonal potential. On this direction many poten- 
tial minima are situated that are occupied by the particles. 
However, an equivalent symmetry direction is found at an 
angle of 27r/5 with the rows. It should also be occupied 
by as much particles as possible. Figure [6] demonstrates 
how this direction can be fitted to the Archimedean-like 
tiling. Colloids close to the red line are marked by red cir- 
cles. In case of a single row of triangles, to move from one 
particle to another close to the red line, one steps Ix/^ 
in x-direction and 2ly in ^/-direction. In case of a double 
row of triangles, the step lenghts are Ix in x-direction and 
3ly in ^/-direction. The sequence of single and double rows 
corresponds to the sequence of the Fibonacci chain. There- 
fore, there are r times as many double rows as single rows 
of triangles, where r = (l + /2^ 1.618 is the number 
of the golden mean. As a consequence, the orientation of 
the red line is determined by 



tan(27r/5) 



(2 + 3r) ly 
{\^r)lx^ 



(5) 



i.e., we find the ratio ly/lx ^ 0.951 which agrees remark- 
ably well with the value determined from simulations in 
subsection 14.11 

To calculate the mean particle spacing a^, we compare 
the number density in the Archimedean-like tiling. 



ay. 



1 



(6) 
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Fig. 7. Average potential energy of a colloid in the decagonal 
potential calculated for a triangular and for an Archimedean- 
like tiling structure as a function of density, which is given by 
the particle spacing as/ay - To determine the average potential 
energy, we constructed a perfect triangular or Archimedean- 
like tiling structure with the corresponding particle spacing 
and calculated the average potential depth at the positions of 
the colloids of the structure. 



with the number density of a triangular lattice p = 2/ (^/Sal) 
By using Ix from Eq. ([5]), we find for the mean particle 
spacing of the Archimedean-like tiling, 



Fig. 8. Average potential energy of a colloid in the 
tetradecagonal potential calculated for a triangular and for 
an Archimedean-like tiling structure as a function of density, 
which is given by the particle spacing as/ay- 



as 



2 2 + 3r 



1 



v/3 ^ + r tan (27r/5) 



_8_ 
15 



(7) 

where ly has to be one of the length scales Ai^2 intro- 
duced in Eq. (|4]). Therefore, Archimedean-like tilings on 
the decagonal substrate occur for 



ds/cLV = 

or as /ay = 



^ 0.579 

0.937. (8) 



Indeed, in experiments and in simulations, we find the 
Archimedean-like tiling structures for densities close to 
these theoretical values (see, e.g., phase diagram in Fig. 

In Fig. [71 we show how the average potential energy 
of a colloid varies with density when it is part of either a 
triangular lattice or an Archimedean-like tiling structure 
under the influence of the decagonal substrate potential. 
We assumed that the Archimedean-like tiling structure is 
stretched perpendicular to the direction of the rows such 
that Eq. (j5j) holds. From Fig.[71we can clearly identify two 
densities where the Archimedean-like tiling fits to the sub- 
strate very well. Most of the colloids are located in minima 
of the potential and therefore the average potential energy 
exhibits two sharp minima. These densities correspond 
to the values calculated in Eq. (|8]). Note, we also deter- 
mined the average potential energy for other Archimedean 
tilings, but did not find any other structure that shows 
pronounced minima. For example, an Archimedean tiling 




Fig. 9. Brownian dynamics simulation without confining 
boundaries started at an initial density that corresponds to 
as /ay = 0.565. (a) Delaunay construction (as in Fig. [J) dis- 
playing the colloidal configuration when the simulation starts 
(t = 0). After time t — 0.05jay /ksT, where 7 is the friction 
coefficient of a colloid, an Archimedean- like tiling structure can 
be observed which remains almost unchanged up to the end of 
our simulation at time t = 0.25jay /ksT shown in (b). Note, 
while the Archimedean-like tilings shown in Fig. [T] are equi- 
librium phases that occur in systems with densities fixed at 
special values, the structure in (b) is only metastable because 
it will slowly dissolve when the density decreases due to the 
missing boundaries. 



of type (3^ • 4 • 3 • 4) has exactly the same particle inter- 
action energy as the tiling of type (3^ -4^). However, by 
comparing the average potential energies, we find that the 
substrate favors the row structure corresponding to type 
(33-42). 

To determine the densities, where Archimedean-like 
tilings occur under the influence of the tetradecagonal sub- 
strate potential, we calculate the average potential energy 
of a colloid within an Archimedean tiling structure as a 
function of the particle spacing (see Fig.[8j). Again, we find 
pronounced minima of the potential energy (although not 
as sharp as for decagonal potentials) at as /ay = 0.557, 
(^s/obv = 0.691, and as /ay = 0.861. At these values we 
were indeed able to observe Archimedean-like tilings in 
the simulations. 
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We just argued that Archimedean- like tilings should 
only occur at specific densities. Correspondingly, when us- 
ing Monte-Carlo simulations to find the equilibrium phases, 
we only observe them close to the calculated densities, 
e.g., in a narrow density region in the phase diagram of 
Fig. [3l Whereas in the simulations the density is a con- 
stant quantity, it can to a certain degree change its value 
in experiments; the experimental system is large enough 
that the density can become non-uniform and colloids can 
leave the boundary box. Indeed, when in the experiment 
we bring the density close to the ones predicted by Eq. ([8j), 
we observe that the colloids self- adjust their density so 
that the Archimedean- like tiling forms. Figure [71 explains 
this observation since the energy is lowered by forming an 
Archimedean-like tiling, at least relative to the triangu- 
lar phase. To further investigate the idea of self-adjusting 
the colloidal density, we performed the following Brow- 
nian dynamics simulations: We confine the colloids at a 
density given by as /ay = 0.565 that is larger than the 
one suitable for Archimedean-like tilings {as /ay = 0.58). 
We then remove the confining boundaries so that the col- 
loids are able to fioat over the substrate potential and the 
system expands. Indeed, we find that particles adjust to 
the potential by forming the characteristic rows of squares 
and triangles of an Archimedean-like tiling (see Fig. [9]). Of 
course, this structure is not the thermodynamic ground 
state in such an open system. However, it remains present 
during the duration of the simulation which is the time 
a free particle would need to diffuse a distance ay- Only 
at the edges of an Archimedean-like tiling domain does 
the structure start to dissolve. This demonstrates while 
Archimedean-like tilings are equilibrium phases at special 
densities only, they also occur as metastable structures in 
systems with open boundaries. 



5 Phasonic Distortions 

In this section we demonstrate, how the ordering of col- 
loids in the Archimedean-like tiling that forms in the decago- 
nal laser field is affected by manipulating the phasonic de- 
grees of freedom as defined in Eq. ([2]). The consequences 
of a constant phasonic displacement on a quasicrystalline 
phase with decagonal symmetry were presented in [19]. 





Fig. 10. Inducing a phasonic flip: Delaunay construction (as in 
Fig. [U for the Archimedean- like tiling phase {as/av = 0.577, 
Vo/{kBT) = 20) obtained by Monte-Carlo simulations for 
Wy = (top) and Wy = O.Olbav (bottom). The flip in the 
Fibonacci chain defined by the single and double rows of tri- 
angles is marked by dotted circles. 



along its five equivalent directions. One can choose pha- 
sonic displacements such that one set of dark parallel lines 
does not change. So if the square and triangular rows are 
parallel to this set of lines, the whole Archimedean-like 
tiling is not affected by the phasonic displacement. For 
example, a phasonic displacement Wx does not infiuence 
the square and triangular rows when they are oriented 
along the x direction. As we show in the next subsec- 
tion, a phasonic drift in Wx may even help to stabilize the 
Archimedean-like tiling structure. 



5.1 Inducing Phasonic Flips 

We first study how a constant phasonic displacement in- 
fluences the Archimedean- like tihng. Figure [TOl shows pat- 
terns of the Archimedean-like tiling in the decagonal laser 
field for different phasonic components Wy = and Wy = 
0.015ay. Note, a whole row of triangles and a row of 
squares have interchanged their position. This corresponds 
to a phasonic fiip in the Fibonacci chain defined by the 
sequence of long and short distances between the square 
rows. 

Figure [2] (d) shows very clearly sets of parallel lines of 
low intensity in the decagonal laser field that are oriented 



5.2 Stabilization by Phasonic Drifts and Gradients 

A steadily growing Wx in time orients the Archimedean- 
like tiling along the x-direction. This is demonstrated in 
Fig. [11] for a phasonic drift velocity Awx = 10~^ay per 
unit time. Domains that are not oriented along the x- 
direction are reduced since a nonzero Awx constantly in- 
duces phasonic fiips in the rows of squares and trian- 
gles of the Archimedean-like tiling unless they are ori- 
ented along the x-direction. Another way of explaining 
the consequences of the phasonic drift velocity Awx goes 
as follows. One can demonstrate that a steadily growing 
phasonic displacement rearranges the decagonal potential 
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Fig. 11. Stabilization by a phasonic drift: Archimedean- like 
tiling phase (same parameters as in Fig. [T0|) obtained by 
Monte-Carlo simulations without phasonic displacement (top) 
and in a potential with a phasonic drift velocity of Awx = 
lO^av per Monte-Carlo step (bottom). The small graph on the 
right-hand side shows the decagonal potential averaged over Wx 
as a function of y. 




Fig. 12. Stabilization by a phasonic gradient: Archimedean- 
like tiling phase (same parameters as in Fig. [T0|) obtained by 
Monte- Carlo simulations without phasonic displacement (top) 
and in a potential with a phasonic gradient with Wx/x = 0.15 
(bottom). 



landscape such that existing potential wells disappear and 
reappear at other locations. Now consider a colloid in the 
decagonal laser field that experiences a sufficiently fast 
phasonic drift in K;aj-direction, i.e., the relaxation time of 
the colloid within a potential well is much larger than 
the time scale on which the potential well vanishes due 
to the phasonic drift. Then one can introduce an effective 
potential {V)w^{y) for the colloids that is the decagonal 
potential V{x^ y) averaged over Wx- It is illustrated at the 
right hand side of Fig. [TTl From Eq. ([1]) it is clear that 
{y)wx{y) does not depend on x or equals V{x^y) aver- 
aged over the spatial coordinate x already introduced in 
Eq. (m, i.e., 

{V)^Ay) = {V)x{y). (9) 

Clearly, such an effective potential can only support do- 
mains of Archimedean-like tilings oriented along the x- 
direction. 

Another possibility of orienting domains of Archimedean- 
like tilings preferentially along the x-direction is to intro- 
duce a phasonic displacement field with a constant gradi- 
ent in x-direction. Figure [13] shows the resulting decagonal 
interference pattern. The gradient in Wx destroys all con- 
tinuous lines of low intensity in the potential landscape, 
except those in x-direction, as indicated by the red lines. 
Therefore, the Archimedean-like tiling is mainly oriented 
along the direction of the phasonic gradient, as illustrated 
in Fig. [H 
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Fig. 13. Interference pattern of the decagonal laser field calcu- 
lated from Eq. ([T]) when the phase 00 attached to wave vector 
Go varies with a constant gradient, 00 = 0.2x/av- Via Eq. 
(|2|) one can show that it creates a constant gradient in the x 
components of the phononic and phasonic displacements. As 
a result, the lines of low intensity are infinitly extended in 
x-direction whereas in all other directions they have a finite 
length (see, e.g., the dark stripes framed by the red lines). 
Similar jags are studied in [24..25j . 
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In the experimental interference pattern of Fig. ^) 
some of the lines of low intensity just end in the middle 
of the pattern. This suggests that phasonic distortions are 
present, as we will now demonstrate. Suppose that the five 
laser beams in Fig. [2^) are not perfectly arranged around 
the vertical. For example, if the tilt angle of beam 
is larger than the other angles, the projected wave vector 
becomes G[) = Go + AGq with \G'q\ > \Gi\ {i = 0, 1...4). 
From Eq. ([T]) one finds that the deviation Z\Go can be in- 
terpreted as a constant gradient in the phase attached 
to Go, (/>o = ^Go • r. Using Eq. ([2]), this phase <po can be 
decomposed into x components of the phononic and pha- 
sonic displacement fields with constant gradients. There- 
fore, laser beams that are not perfectly adjusted relative 
to the vertical lead to a substrate potential with phononic 
and phasonic distortions. Furthermore, such distortions 
can also be controlled intentionally by displacing one of 
the laser beams radially from its position on an ideal pen- 
tagon. This way one can create a preferred orientation for 
the Archimedean-like tiling, as discussed above. 

In summary, phasonic drifts or distortions can be em- 
ployed to grow large domains of an Archimedean-like tiling 
structure with a preferred orientation. These domains are 
stable when the phasonic drift is stopped or the phasonic 
distortion is removed. 



6 Outlook and Conclusions 

In this article we have studied the structure of Archimedean- 
like tiling phases that form in decagonal and tetradecago- 
nal laser fields at certain densities. We have also demon- 
strated how the domain size of such tilings is influenced by 
phasonic distortions. In particular, if the interfering laser 
beams are not adjusted properly, phasonic displacement 
fields with a constant gradient arise. Interestingly, simu- 
lations using substrate potentials with further quasicrys- 
talline symmetries also show similar phases. For example, 
in an interference pattern of eight beams we find rows 
of squares and triple or four-fold rows of triangles [Fig. 
[TUf a)]. A laser field created by ten beams favors exactly 
the same structure as the decagonal potential obtained by 
five laser beams [Fig.fMlfb)]. For the case of eleven beams, 
we also observe the characteristic rows of triangles and 
squares [Fig. [TWc)]. However, we are not able to analyze 
this ordering in more detail with respect to the order per- 
pendicular to the rows, since the domains are too small. 
It is an interesting question for future research to under- 
stand, why such Archimedean-like tilings form at specific 
densities on substrates with different quasicrystalline sym- 
metries. Depending on the specific rotational symmetry, 
all these structures consist of single and double (or even 
triple or four-fold) rows of triangles or squares that ulti- 
mately form a quasiperiodic order in one direction. 

While on the decagonal substrate the rows in the Ar- 
chimedean-like tiling form a Fibonacci chain, they show 
domains with periodic order on the tetradecagonal sub- 
strate interrupted by double rows of squares. We find that 
these double rows, as the double rows of triangles on the 




Fig. 14. Delaunay construction (as in Fig. [T]) of simulation 
results that consist of rows of triangles and rows of squares in 
laser fields (of strength Vo/(/csT) = 50) created by (a) eight 
beams {as /ay = 0.8), (b) ten beams {as /ay = 0.7), and (c) 
eleven beams {as /ay — 0.74). 

decagonal substrate, occur at symmetry centers of the sub- 
strate potential. As discussed in [20^, the number of the 
symmetry centers is by a factor of 100 larger in the decago- 
nal compared to the tetradecagonal substrate. Therefore, 
they enforce the formation of a Fibonacci chain in the first 
case, whereas in the tetradecagonal laser field larger pe- 
riodic domains form due to the much smaller number of 
symmetry centers. 

Amongst other phases (see Ref. [19 ), Archimedean- 
like tilings appear as a compromise between structures 
favored either by the particle interaction or by the sub- 
strate potential. There is some evidence that they occur 
in experiments that study how atoms order on atomic sur- 
faces with decagonal symmetry [18] . For the specific exam- 
ple of the Archimedean-like tiling, we have presented here 
Brownian dynamics simulations to demonstrate how the 
colloidal adsorbate reacts on phasonic distortions in the 
substrate. This rises the question if phasonic excitations 
are detectable in atomic substrates. The Archimedean-like 
tiling as an adsorbate phase might be a good candidate 
to visualize these excitations in the underlying substrate 
since resulting rearrangements in the rows of triangles and 
squares should be visible. 

We would like to thank R. Lifshitz and H.-R. Trebin for help- 
ful discussions. We acknowledge financial support from the 
Deutsche Forschungsgemeinschaft under Grant No. RO 924/5- 
1 and BE 1788/5. 
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